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Abstract 

In this paper, we consider a general bi-directional relay network with two sources and N relays 
when neither the source nodes nor the relays know the channel state information (CSI). A joint relay 
selection and analog network coding using differential modulation (RS-ANC-DM) is proposed. In 
the proposed scheme, the two sources employ differential modulations and transmit the differential 
modulated symbols to all relays at the same time. The signals received at the relay is a superposition 
of two transmitted symbols, which we call the analog network coded symbols. Then a single relay 
which has minimum sum SER is selected out of N relays to forward the ANC signals to both 
sources. To facilitate the selection process, in this paper we also propose a simple sub-optimal 
Min-Max criterion for relay selection, where a single relay which minimizes the maximum SER 
of two source nodes is selected. Simulation results show that the proposed Min-Max selection has 
almost the same performance as the optimal selection, but is much simpler. The performance of the 
proposed RS-ANC-DM scheme is analyzed, and a simple asymptotic SER expression is derived. 
The analytical results are verified through simulations. 
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I. Introduction 

In a bi-directional relay network, two source nodes exchange their messages through the 
aid of one or multiple relays. The transmission in bi-directional relay network can take place 
over either four, three or two time slots. In the four time slots transmission strategy, the relay 
helps to forward source Si's message to source 5*2 in the first two time slots and source 
S'2's message to source Si in the next two time slots. Four time slots transmission has been 
shown to be very inefficient. When the relay receives two sources' messages, it combines 
them before forwarding to the destination, which will save one time slot transmission. This 
three time slots transmission scheme is usually referred to as the digital network coding [IJ- 
Q. In this method, two source nodes transmit to the relay, separately. The relay decodes the 
received signals, performs binary network coding, and then broadcasts it back to both source 
nodes. 

To further improve the spectral efficiency, the message exchange between two source nodes 
can actually take place in two time slots. In the first time slot, both source nodes transmit at 
the same time so that the relay receives a superimposed signal. The relay then amplifies the 
received signal and broadcasts it to both source nodes in the second time slots. This scheme 
is referred to as the analog network coding (ANC) ||4l-[[6l. Various transmission schemes 
and wireless network coding schemes in bi-directional relay networks have been analyzed 
and compared in [TTll- lfTljl . 

Most of existing works in bi-directional relay communications consider the coherent detec- 
tion at the destination and assume that perfect channel state information (CSI) are available 
at the sources and relays IfTTl- lfTll . In some scenarios, e.g. the fast fading environment, the 
acquisition of accurate CSI may become difficult. In this case, the non-coherent or differential 
modulation would be a practical solution. In a differential bi-directional relay network, each 
source receives a superposition of differentially encoded signals from the other source, and 
it has no knowledge of CSI of both channels. All these problems present a great challenge 
for designing differential modulation schemes in two-way relay channels. 

To solve this problem, in [[T3l . a non-coherent receiver for two-way relaying was proposed 
for ANC based bi-directional relay networks. However, the schemes result in more than 3 
dB performance loss compared to the coherent detection. To further improve the system 
performance, a differential ANC scheme was proposed in [[T4l and a simple linear detector 
was developed to recover the transmitted signals at two source nodes. Results have shown 
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that it only has about 3 dB performance loss compared to its coherent counterpart. 

Recently, it has been shown that the performance of wireless relay networks can be further 
enhanced by properly selecting the relays for transmission ifTSl - lfTTl . Consequently, it is 
beneficial to design an effective relay selection scheme for the bi-directional transmission 
scheme with multiple relays as well in order to achieve spatial diversity. In this paper, we 
propose a bi-directional joint relay selection and analog network coding using differential 
modulation (RS-ANC-DM) so that the CSI is not required at both sources and the selected 
relay. In the proposed BRS-DS-ANC scheme, two source nodes first differentially encoder 
their messages and then broadcast them to all the relays at the same time. The signals 
received at the relay is a superposition of two transmitted symbols. Then a single relay 
which minimizes the sum SER of two source nodes is selected out of all relays to forward 
the ANC signals to both sources. Each source node then performs the differential detection 
and subtract its own message to recover the message transmitted by the other source node. 

The performance of optimal relay selection is very difficult to analyze. To facilitate the 
analysis and selection procedure, in this paper we also propose a simple sub-optimal Min-Max 
criterion for relay selection, where a single relay which minimizes the maximum BER of two 
source nodes is selected. Simulation results show that the proposed Min-Max selection has 
almost the same performance as the optimal selection, but is much simpler. The performance 
of the proposed BRS-ANC scheme is analyzed, and an asymptotic SER expression is derived. 
The analytical results are verified through simulations. 

The rest of the paper is organized as follows: In Section II, we describe the proposed BRS- 
DS-ANC scheme. The performance of the RS-ANC-DM is analyzed Section III. Simulation 
results are provided in Section IV. In Section V, we draw the conclusions. 

Notation: Boldface lower-case letters denote vectors, (■)*, (■)^ and (■)^ represent con- 
jugate, transpose, and conjugate transpose, respectively. E is used for expectation, Var 
represents variance, ||xp = x^x, and denotes real part. 

II. Joint Relay Selection and Analog Network Coding using Differential 

Modulation 

We consider a general bi-directional relay network, consisting of two source nodes, denoted 
by Si and S2, and relay nodes, denoted by . . . ,Rn- We assume that all nodes are 
equipped with single antenna. In the proposed RS-ANC-DM scheme, each message exchange 
between two source nodes takes place in two phases, as shown in Fig. [TJ In the first phase. 
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both source nodes simultaneously send the differentially encoded information to all relays 
and the signal received at each relay is a superimposed signal. In the second phase, an optimal 
relay node is selected to forward the received signals to two source nodes and all other relay 
nodes keep idle. 

A. Differential Encoding and Decoding in Two-Way Relay Channels 

Let Cj(t), i = 1,2, denote the symbol to be transmitted by the source Si at the time t. We 
consider a MPSK modulation and assume that Cj(t) is chosen from a MPSK constellation of 
unity power A. Source i first differentially encodes the information symbols Ci{t) 

Si{t) = Si{t - l)Ci{t) (1) 

The differential encoded signals are then simultaneously transmitted by two source nodes 
with unit transmission power to all the relays. The signal received in the k-th relay at time 
t can be expressed as 

yr,k{t) = hl,kSl(t) + h2,kS2{t) + nr,k{t), (2) 

where hi^k, ^ = 1, 2, = 1, A^, is the fading coefficient between Si and Rk- hi this paper, 
we consider a quasi- static fading channel for which the channels are constant within one 
frame, but change independently from one frame to another, n^^kit) is a zero mean complex 
Gaussian random variable with two sided power spectral density of Nq/2 per dimension. 

Upon receiving the signals, the relay Rk then processes the received signal and then 
forwards to two source nodes. Let Xr^k{t) be the signal generated by the relay R^ and it 
is given by 

Xr,k{t) = I3ky*.^kii)^ (3) 

where (3k is an amplification factor, so that the signal transmitted by the relay satisfy the 
following power constraint 

E{\Xr,kit)\^) = 1. (4) 

We should note that unlike the traditional ANC schemes H-im, the relay forwards the 
conjugate of the received signal. The reason of doing this is to facilitate the differential 
detection at the destination fT5|. 

Substituting Eqs. ^ and ([3]) into ([H), we can derive (3k, 



^'~''\h.k\' + \h.k\' + No 
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However, as the relay has no CSI, /3fc has to be obtained in other way. Let y^;. = 

[?/,.fc(l),---,2/r,fc(^)]^, Si = [Si(l),...,Si(L)]'^, S2 = [S2(1),...,S2(L)]^, n,.fe = K,fc(l), • • • 

where L is the frame length. Then we can rewrite the received signals in ([2l) in a vector format 
as follows 

yr,k = hl,kSl + /l2,fcS2 + n^,fc, (6) 

and (3k can be then approximated by the k-th relay node as 



Pk = \l ^ — ^V^' 

After deriving (3k, the relay Rk then forwards Xr^it) to two source nodes. Since Si and 5*2 
are mathematically symmetrical, for simplicity, in the next we only discuss the decoding as 
well as the analysis for signals received by ^i. The signal received by at time t, denoted 
by Vi^kit), can be written as 

= (3khi,kyr,ki't) +ni,kit) 

= fJ'kSlit) + UkSlit) + Wi^kit) 

= iJkslit) + Uksl{t - l)cl{t) + wi^k{t), (8) 

where ^k = (3k\hi±\^ > 0, Uk = (3khi^kh*2k^ and wi,fe(t) = (3khi^kn*k(t) + rii^kit). 

Since Si{t) is known to the source 5*1, to decode C2(t), Si needs to estimate fik and z/^. 
Let y^ k = • • • and Wi,^ = . . .,wi,kiL)f- Then at high SNR, we 

can obtain the following approximation 

f^l + Wll-^^-^- (9) 
Since the source node Si can retrieve its own information si{t — 1) and ci{t), we have 

yi,ki^) - c\{t)yi^kit - 1) - yi,kit) 

= VkSl{t - 1) (Ci(t) - C2(t))* + Wi,k{t), (10) 

where Wi^kit) — Ci{t)wi^kit — 1) + Wi,fc(^)- Then, |z/fcp can be approximated as [fT4l 

~H ~ 

I |2 ^ yi,fcyi,fc (\\\ 

i'''^' ^ m\\ci{t)-cm^X ^''^ 

where y^^k = [yi,A:(l), • • • ,yi,kiL — 1)]^, and E[|ci(t) — C2(t)p] is a constant which can be 
pre-calculated by two source nodes, which is given in Appendix A. As fik is positive, it can 
be derived from ^ and (fTTl) as 
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/^fc ~ (e,)+ (12) 

where 6^ = — |z/jtp, (X)+ is equal to X when X > and otherwise is equal to 0. 

By subtracting /ifcS^(t), ([8]) can be further written as 

= VkSl{t - + Wl,fc(t) 

= {y[^^{t - 1) - - 1)) cl{t) + (13) 

Finally, the following linear detector can be applied to recover C2(t) 

C2{t) = arg max Re {y'l^{t - l)y;,,(t)c2(t)} . (14) 



B. Relay Selection with Differential Modulation in Two-Way Relay Channels 

In the proposed RS-ANC-DM scheme, only one best relay is selected out of relays to 
forward the received ANC signals in the second phase transmission. We assume that at the 
beginning of each transmision, some pilot symbols are transmitted by two source nodes to 
assist in the relay selection. One source node (either source Si or 5*2) will determine the one 
best relay according to a certain criterion and broadcast the index of the selected relay to 
all relays. Then only the selected relay is active in the second phase of transmission and the 
rest of relays will keep idle. 

1 ) Optimal Single Relay Selection: For the optimal single relay selection, a single relay 
node, which minimizes the sum SERs of two source nodes, i.e. SERi ^ + SER2,a:, will be 
selected, where SERi ^ and SER2,fc represent the SERs at source nodes and 5*2, respectively. 
The main challenges in relay selection for differential modulation as mentioned before is that 
the relay node is determined by only one source (either Si or S2) without the knowledge of 
any CSI. For simplicity and without loss of generality, in the next we use Si to select the 
optimal relay node. Obviously, the main difficulty here is to estimate SER2,fc. 

For M-PSK constellations, the conditional SER assoicated with the A;-th relay at the source 
Si is given by [UHl 

(AJ~l)7r 

SERi,fc(/ii,,, h2,k) = - [ " exp f 3ps^.l'i,,k \ ^^3^ 

where 7^1 ^ is the effective SNR at the source Si and gp^^ = sin^ jj. As CSI is unknown to 
the receiver, the effective SNR ^di,k has to be estimated without knowledge of CSI. 
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By ignoring the second order term, the corresponding SNR of the proposed differential 
detection scheme in (fT3]) can be written as 



Var{2w;i,fc(t)} 





|2 
1,A:| 






\il'2,k\ 




? +^s\h2,k\'^ 



(16) 

where Var{u;i,fe(t)} ^ /3^A/'o|/ii,A:P + A^o, V's = 4^, and = 

Recalling /ife = Vk — f^khi^kh^^, and their corresponding estimates in (fTTl) and 

(fT2l) . 7rfi,fc in (fT6l) can be further calculated as 

1^ |4 1 ^ 1 2 

2(2|/Xfc|2+|z/,|2)(K|2+|z/,|2)iVo' ^^^^ 

Similar to (fTTl) . the SNR of the proposed differential detection scheme in the source 5*2 
can be written as 

ld2,k 









2,fc 


2 


2 


A\h2,k 




^k\ 


2(2|z/fc| 




P+kA:|2)iVo 



(18) 

And its SER can be then calculated as 

(M-l)w 

SER2,,(/ii,.,M = - / " exp(-^^^^)de. (19) 
7T Jq V sm / 

Among all relays, the destination will select one relay, denoted by TZ, which has the 

minimum destination SER: 

7^ = min {SERi,fc(/ii,fc, h2,k) + SER2,fe(/ii,fc, h2,k)} ,kel,...,N (20) 

k 

2) Sub-Optimal Single Relay Selection: The optimal single relay selection scheme de- 
scribed in the above section is very difficult to analyze. In this section we propose a sub- 
optimal single relay selection scheme. It is well-known that the sum SERs of two source 
nodes (SERi ^ + SER2,a,.) is typically dominated by the SER of the worst user. As a result, for 
low complexity, the relay node, which minimizes the maximum SER of two users, can be 
selected to achieve the near-optimal SER performance. We refer to such a selection criterion 
as the Min-Max selection criterion. Let TZ denote the selected relay. Then the Min-Max 
selection can be formulated as follows, 

TZ = minmax{SERi,fc(/ii,fc, h2,k),SER2,k{hi,k, h2,k)} ,k e 1, . . . ,N, (21) 

k 
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which can be further formulated by using the effective SNRs 

7^ = maxmm{'ydi,k,ld2,k},k el,...,N, (22) 

k 

where the calculation of ^ and 'jd2,k can be obtained from (flTI) and (fTSi) . respectively. 

III. Performance Analysis of RS-ANC-DM scheme Based on Min-Max 

Selection Criterion 

In this section, we derive the analytical average SER of the proposed differential bi- 
directional relay selection schemes. As mentioned before, the optimal relay selection scheme 
is very difficult to analyze. As it will be shown later, the Min-Max selection scheme proposed 
in section II has almost the same performance as the optimal selection scheme. Therefore, 
in this section, we will analyze the RS-ANC-DM scheme based on the Min-Max selection 
criterion. 

A. Sub-Optimal Relay Selection 

For the Min-Max selection criterion, the effective SNR of the selected relay TZ can be 
expressed as follows, 

77e = maxmin{7d^,fc, 7d2,fc}, k e 1,...,N. (23) 

Now let us first calculate the PDF of 77^. As 'jdi.k and 7^2 .a: are identically distributed, 
they have the same PDF and CDF, denoted by f'y^.{x) and Fj^{x), respectively. 

Define Yk^^—''^^^{'ydi,k,'yd2,k}- Let /^min(x) and F^mi.i(x) represent its PDF and CDF, 
respectively. For simplicity, assuming that 7^^ ^ and 7^3 ^ are independent, then the PDF of 
TZ can be calculated by using order statistics as |fT9l 

f,Jx) = Nf^^.nix)F^p{x) = 2Nf,,ix)il - F,,ix))[l - (1 - F,,ix))Y-\ (24) 

where f^^in{x) = 2/^^(x)(l — F^^{x)), F^mm(x) = 1 — (1 — F^^{x)Y, and f^^{x) can be 
found in IEQII 



2xexp {—x{ip^ + )) 



(25) 

where Kq{-) and Ki{-) denote the zeroth-order and first-order modified Bessel functions 
of the second kind, respectively, and U{-) is the unit step function. At high SNR, when 
z approaches zeros, the Ki{z) function converges to 1/z [21], and the value of the Kq{-) 



U{x) 
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function is comparatively small, which could be ignored for asymptotic analysis. Hence, at 
high SNR, fj^^ix) in (l25l) can be reduced as 



/7fc(^) ~ ^exp f-^a; ) , (26) 



where ip = 2{ijj^ ^ + ip^^). Its corresponding CDF can be written as 



2 

The PDF of 77^ can then be approximately calculated as 



F^,(a;)«l-exp(-^x). (27) 



/^^ (x) ^ Nip exp {-ipx) [1 - exp{-ipx)f~^ . (28) 

The CDF of 77^ can be approximately written as 

F^^{x) ^ [1 - exp{-iljx)f ^ {tjjxf , (29) 

where lim 1 — exp(— x) = x. 

x—^oo 

The SER conditioned on the instantaneous received SNR is approximately [|22l 

SER(77e|/ii,/i2)~g(v^), (30) 

where Q{-) is the Gaussian-Q function, Q{x) = exp(— t^/2)dt, and c is a constant 

determined by the modulation format, e.g. c = 2 for BPSK constellation. 

The average SER can be then derived by averaging over the Rayleigh fading channels by 
SER(77^) = E [SER{^n\h, /12)] =^[Q (v^)] • (31) 

By introducing a new random variable (RV) with standard Normal distribution X ~ 
A/^(0, 1), the average SER can be rewritten as [|23l 

SER(77e) = P{X>y^} 

r x2 

= P <^ 77^ < — 



E 



X 



c 



^\ f 2N f -^^ 



27r V c / in V 2 



x^'^exp -— dx. (32) 



Based on the fact that f^^ t'^'' exp{-kt'^)dt = fey^v^ [El, we can finally obtain 



2(2fc)" 

(2iV-l)!! ftp 



N 



SER(77e) = ^ r^ic' ' ^^^^ 



(2n-l)! 
n!2" ■ 



where (2n - 1)!! = HLi 2^ - 1 

It clearly indicates in (l33l) that a diversity order of can be achieved for the proposed 
RS-ANC-DM scheme in a bi-directional relay network with two sources and N relays. 
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IV. Simulation Results 

In this section, we provide simulation results for the proposed RS-ANC-DM scheme. We 
also include the corresponding coherent detection results for comparison. All simulations 
are performed with a BPSK modulation over the Rayleigh fading channels, and the frame 
length is L = 100. In order to calculate Uk in (fTTl) . we use the normal constellation and the 
constellation rotation approaches introduced in Appendix A. For simplicity, we assume that 
Si, S2, and Rk (k = 1, . . . , N) have the same noise variance. 

A. Simulated Results 

From Fig. [2] to Fig. HI we present the simulated SER performance for the proposed RS- 
ANC-DM schemes. The performance of the corresponding coherent detection are plotted as 
well for better comparison. In Fig. [21 we compare the optimal relay selection method and 
the sub-optimal Min-Max relay selection method. It can be observed from the figure that 
the proposed Min-Max selection approach has almost the same SER as the optimal one. In 
particular, when the number of relay nodes increases, we almost cannot observe any difference 
between these two methods, which indicates that the Min-Max relay selection achieves near 
optimal single relay selection performance. 

In Fig. [3l we compare the RS-ANC-DM scheme with its coherent detection counterpart. 
It can be noted that the differential scheme suffers about 3 dB performance loss compared to 
the corresponding coherent scheme. We can also see from Fig. |3] that the SER performance 
is significantly improved when the number of the relay increases. In Fig. H we include the 
Genie-aided result by assuming that n is perfectly known by the source such that traditional 
differential decoding can be performed. It shows from the results that there is almost no 
performance loss using the estimation method in (fT2l) which clearly justifies the robustness 
of the proposed differential decoder. 

Fig. |5] compares the simulated SER performance for our proposed RS-ANC-DM and 
the non-coherent schemes [fT3]| in bi-directional relaying without using constellation rota- 
tion, where N = 1,2,4,8. It can be observed that our proposed scheme has much better 
performance than the detector in [fT3l . The main reason is that the non-coherent detection 
approach employed in [[T3]| statistically averages off the impact of channel fading coefficients 
by ignoring the instantaneous channel state information and thus causes much performance 
loss. Comparatively, our proposed differential detection is a symbol by symbol based detection 
and is thus be able to adapt to the variation of the channel. 
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B. Analytical Results 

In Fig. [6l we compare the analytical and simulated SER performance of the proposed 
differential modulation scheme. I n order to obtain fine estimation in (fTTI) . the signal constel- 
lation used by Si is rotated by 7r/2 relative to that by 5*2. From the figure, it can be observed 
that at high SNR, the analytical SER derived by (|33l) is converged to the simulated result 
using optimal relay selection. This verifies the derived analytical expressions. 

C. Constellation Rotation 

In Fig Ul we examine the SER results of the proposed differential modulation scheme 
in comparison with the one without using constellation rotation, as shown in Appendix A, 
where the signal constellation used by 5*1 is rotated by n/2 relative to that by 5*2. It can be 
observed that the new result has very similar with the curve without rotating constellations. 
This indicates that using constellation rotation may not obtain any gains given large frame 
length. 

V. Conclusions 

In this paper, we proposed a joint relay selection and analog network coding using differ- 
ential modulation over two-way relay channels when neither sources nor the relay has access 
to the channel state information. A simple Min-Max relay selection method is proposed and 
it has been shown that it achieves almost the same performance as the optimal single relay 
selection scheme. An asymptotic SER expression is derived. It is shown that the proposed 
RS-ANC-DM scheme can achieve the full diversity order of for the system with relays. 
Results are verified through simulations. 

Appendix A 
The Calculation of E[|ci(t) - C2(t)p] in (fTTI) 

From (fTTI) . it shows that the average power of Ci(t) — C2(t) needs to be calculated. When 
M-PSK constellations are applied, the number of symbols produced in the new constellations 
by ci(t) — C2(t) is finite. Hence, it is easy to derive the average power of the new constellation 
sets. We refer to this as the normal constellation approach. 

Note that the value of ci(t) — C2{t) can be equal to zero, which may affect the estimation 
accuracy in (fTTI) . In order to overcome this problem, we may properly choose a rotation 
angle for the symbol modulated in source ^2 by C2{t)e~^^, ensuring that Ci(t) — C2{t) in 
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(fTTI) is nonzero. For a M-PSK constellation, the effective rotation angle is in the interval 
[— 7r/Af, tt/M] from the symmetry of symbols. For a regular and symmetrical constellation, 
the rotation angle may be simply set as = n /M. Similar approach may be used to generate 
the rotation angle for other types of constellations. 
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Phase 1 : Transmission via orthogonal channels 




Phase 2: Broadcasting via orthogonal channels 



Fig. 1. Block diagram of the proposed BRS-DANC scheme. 




Fig. 2. Simulated SER performance by optimal and sub-optimal detections. 
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Fig. 3. Simulated SER performance by differential and coherent detections, where = 1, 2, 4, 8. 




Fig. 4. Simulated SER performance by differential and Genie-aided detections, where TV = 1, 2, 4, 8. 
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Fig. 5. Simulated SER performance comparisons by our proposed differential approach and the non-coherent scheme 
in [m, where iV = 1,2,4, 8. 
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Fig. 6. Analytical and Simulated SER performance by the proposed differential scheme, where A*" = 1, 2, 4, 8. 
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Fig. 7. Simulated SER performance by the proposed differential scheme with and without using constellation rotation, 
where = 1,2,4,8. 



